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We address a problem of identifying a given pure state with one of two reference pure states, 
when no classical knowledge on the reference states is given, but a certain number of copies of them 
are available. We assume the input state is guaranteed to be either one of the two reference states. 
This problem, which we call quantum pure state identification, is a natural generalization of the 
standard state discrimination problem. The two reference states are assumed to be independently 
distributed in a unitary invariant way in the whole state space. We give a complete solution for 
the averaged maximal success probability of this problem for an arbitrary number of copies of the 
reference states in general dimension. It is explicitly shown that the obtained mean identification 
probability approaches the mean discrimination probability as the number of the reference copies 
goes to infinity. 

I. INTRODUCTION 

Suppose we are given an unknown quantum pure state p on a d dimensional vector space . We know that the 
input state p is either one of two reference states pi and p2, each also being a pure state on C^. What is the best 
strategy to identify the input state with one of the two reference states? 

We can consider two cases depending on what kind of information on the reference states is available. In the first 
case, we are given complete classical knowledge on the reference states pi and p2- This is the problem of quantum 
state discrimination, which was solved by Helstrom 

On the other hand, we can also consider the case where only a certain number {N) of copies of pi and p2 arc given, 
with no classical knowledge on them available. In this case, we could obtain only limited classical information on the 
reference states, since the no-cloning theorem does not allow us to increase the number of copies of the reference 
states. The best we can do is that we perform a POVM measurement on the total state p ^ pf^ (g) pf^ and try to 
identify the input state p with one of the reference states pi and p2. In this paper, this problem is called "quantum 
state identification" . If the number of copies N is infinite, the problem is reduced to quantum state discrimination. 
This is because we can always obtain complete classical knowledge of a quantum state if we have infinitely many 
copies of the state. Thus quantum state identification is a natural generalization of quantum state discrimination, 
which is one of the most fundamental problems in quantum information theory. 

In the case of qubit (d = 2), similar problems but in different setups have already been studied 0,0, 13- Sasaki et 
al. studied quantum matching problem i, where a certain number oi copies of input and reference states are given 
according to a known a priori independent probability. Using the fidelity as a figure of merit, they determined the 
optimal fidelity in the case where the input and reference qubits are independently distributed on a great circle of the 
Bloch sphere (see also Ref. Q)- Bergou et al. Q recently discussed a generalization of unambiguous discrimination 
problem In this generalization we are given three qubit systems, where one of the three qubits, input system, 

is guaranteed to be prepared in one of the states of the other two reference qubits with some occurrence probability. 
Without any classical knowledge on the reference states, we are to tell unambiguously in which of these reference states 
the input qubit is prepared. In Ref. Q the optimal success probability was derived as a function of the occurrence 
probability of the two cases. 

In this paper, we assume that the input state p is guaranteed to be prepared in one of the two reference states pi 
and p2 like in the original discrimination problems and the generalization of unambiguous discrimination discussed in 
Ref. . The two reference states are assumed to be independently distributed on the whole d dimensional pure state 
space C"^ in a unitary invariant way, which will be precisely defined in the next section. No classical knowledge on 
the reference states is available, but only a certain number (N) of their copies are given. Our task is to successfully 
identify the input state with one of the reference states. 

We will give a complete solution of optimal strategies and the mean success probability as a function of the number 
of copies N of the reference states and the dimension d of the state space. We also study the large N limit of the 
mean identification probability and verify that it approaches the mean discrimination probability as the number of 
copies N goes to infinity. 
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II. AVERAGED MAXIMUM PROBABILITY OF PURE STATE DISCRIMINATION 

First we calculate the averaged maximum discrimination probability in order to compare it with the maximum 
identification probability derived later. The success probability of discrimination between two known pure states pi 
and p2 is given by 



1 ^ 

P{PliP2) = l^^tr[EaPa\, (1) 



a=l 



where the occurrence probabilities of pi and p2 are assumed to be equal. Here E2] is a set of POVM and with 
outcome labeled by a{= 1, 2) the input state is guessed to be pa- The POVM may depend on the reference states pi 
and y02, since we have the classical knowledge on pa in this discrimination problem. 

For a given set of pure states pi and p2 , the maximum success probability is expressed in terms of the trace distance 
between the two states as 

Pmax(pi,P2) = ^(1 + £'(P1,P2)), (2) 

where D{pi,p2) = |tr|pi — P2I Q- We note that the trace distance between two pure states is given by 



D{pi,p2) = Vl-tr[PiP2]- (3) 

Now we specify the distribution by which the reference states pi and p2 are chosen. We assume that the two states 
are chosen independently and uniformly from the pure state space on in a unitary invariant way. Denoting the 
average over this distribution by <•••>, we define the averaged maximal success probability of state discrimination 

to be Pmax(d) =< Pmax(/5l,P2) >• 

More precisely, the average is defined as follows. Let us expand a normalized pure state | G by an orthonormal 
reference base {|i)}f=i Sts \(f>) = X^i'^*!*) with J^i = 1- The distribution of \(p) is assumed to be given by the 
uniform distribution of {xi = SRci, yi = '^Ci}f^i on the 2d — 1 dimensional hyper sphere. 

When averaging the distance D{pi,p2), we can assume the state pi is given by 1 1, where | 1 ) is one of the 
vectors in the reference base, which is justified by the unitary invariance of the distance D and the distribution of 
states. 

Expanding the state p2 in the reference base with coefficients c;, we can express the averaged distance as 



fdcdc+Jl^\ci\^ 



where the integration measure corresponding to the uniform distribution on the hyper sphere is given by 

d 



dcdc+ EE / llidx,dy,)5iJ2ix^ + yf) - 1). (5) 

i=i i 



Performing the above integrations, we find that the average distance takes the form 



<D{pi,p2)>^^—^, (6) 
d- 2 

which gives the averaged maximum success probability of discrimination to be 

Prn..{d) = - + (7) 

We note that the mean distance < D(pi,p2) > and consequently the mean identification probability Pinax(rf) 
approach unity as the dimension d goes to infinity. 



III. QUANTUM STATE IDENTIFICATION: CASE OF iV = 1 



Let us consider the simple case, N = 1, where only one copy of each reference state is given. We have three systems 
numbered by 0,1, and 2, each on C"^. The input state p is assumed to be given in system and the reference states 
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pi and p2 in system 1 and system 2, respectively. We specify the system that an operator acts on by a number 
in the parentheses; p{Q) is a density operator on system for example. Our task is then to distinguish two states, 
/9i(0)pi(l)/92(2) and p2(0)/3i(l)yO2(2) , which are assumed to occur with equal probabilities. The POVM i?2} is a 
set of operators on the total system (g) 1 (g) 2 and should be independent of pi and p2 , since no classical knowledge 
on pi or p2 is available in the state identification problem considered in this paper. 
The averaged success probability of state identification is then expressed as 

P^^^'Hrf) = ^ E (t^- {EaPaiO)piil)p2{2)] ), (8) 

a=l,2 

where the average < • • • > should be performed by the unitary invariant and independent distribution of the two 
reference states with the integration measure defined in Eq.(|Sl. 

Here it is very helpful to use the following formula for the unitary average of the tensor product of n identically 
prepared pure states 0: 

< P^" f^, (9) 

where 5„ is the projector onto the totally symmetric subspace of {C"^}®" and the dimension of the subspace is given by 
dn = tr [Sn] ~ n+d~\C d-\- This relation can be derived either by an explicit computation with the integration measure 
(O or more simply by the following group theoretical argument. For any unitary \J on C', we have = 
(p®") , implying that the operator (p®") on the totally symmetric subspace of {C"^}®" commutes with /7®" for any 
V . Shur's lemma then requires that (p®") be proportional to 5„, since J/®" acts on the totally symmetric space 
irreducibly. The proportional coefficient turns out to be l/c?n by a trace argument. Thus we obtain formula 

Performing the average integration by the use of formula Q and using the completeness of POVM, E2 = 1, 

we find 

P^''='Hd) = 1 + ^tr [iJi (52(01) - 52(02))] , (10) 

where ^2(01) and 52(02) are the projectors onto the totally symmetric subspace of space 0® 1 and 0® 2, respectively. 
For any Hermitian traceless operator D, we observe the relation 

max tr\ED] = itrlZ?!. (11) 
o<£;<i ^ ^ 2 ' ' ^ ^ 

Here the maximum is attained if and only if i? = Vd+ + Eu^, where Vd^ is the projector onto the eigenspace 
with positive eigenvalues of D, and Ed^ is any operator in the subspace of zero-eigenvalue of D with the condition 

< Edo < 1- Note that the operator E can be chosen as a projector by taking E^g = Vo, where Vo is is a projector 
onto a subspace spanned by some zero-eigenvalue states of D, which leads to a projective optimal measurement. 

Taking D = 52(01) - 52(02), we obtain 

pil^'Hd) = 1 + ^tr|52(01) - 52(02)1. (12) 

The optimal POVM can be given by the projective measurement, Ei = Vd+ + Vq and E2 = 1 — Ei. 

In order to evaluate p^ax^\d) further, we must determine non-zero eigenvalues and their multiplicities of the 
operator D = 52(01) — 52(02) = i(T(01) — r(02)), where T(01) is the transposition operator between systems and 

1 and T(02) is defined similarly. 

It is instructive to explicitly work out the case of d = 2 first. It is easy to see that the totally symmetric states are 
eigenstates of D with zero eigenvalue 

£)|000) = 0, 
D|lll) =0, 

i:)(|100) + |010) + |001)) = 0, 

i:»(|oii) + |ioi) + |iio)) = o. (13) 

For mixed symmetric states, we find four eigenvectors in the following way. 

D\±: 100) = ■ 100), 

i^|±:110) = ±^|±:110), (14) 
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where 

|±:100) ^ (-2|100) + (1tV3)|010) + (1±V3)|001))/VT2, 

|±:110) = (-2|011) + (lTV3)|101) + (l±V3)|110)j/VT2. (15) 

There is no totally antisymmetric state in this d — 2 case. Thus eigenvalues of D are four O's, two ^'s and two — -^'s, 
which yields tr|I?| — 2^/3. The optimal measurement in this case is given by the following projective measurement: 

El = 1+ : 100)(+ : 100| + 1+ : 110)(+ : 110| +7^0, 

E2 = I- : 100)(- : lOOl + l- : 110>(- : 1101+53(012)- -Po, (16) 

where 53(012) is the projector onto the totally symmetric space and Vq is a projector onto any subspace of it. One 
might wonder why the freedom in choosing Vq does not affect the optimal success probability. The reason is that if the 
whole system is measured and found to be totally symmetric, we conclude that the system was prepared equally likely 
in one of the two states, pi(0)pi(l)p2(2) or p2(0)pi(l)p2(2), and any guess does not improve the success probability. 

Rather surprisingly, this property of the spectrum of operator D in the qubit case prevails in general dimensions d. 
This can be most conveniently seen by computing Z?^ as follows: 

= ^ (2 - T(01)r(02) - r(02)T(01)) = ^ (1 - 53(012) - ^3(012)) = ^M3(012), (17) 

where ^3(012) and A^3(012) are the projectors onto the totally antisymmetric and the mixed symmetric subspace of 
space 18) 1 (8) 2, respectively. From this we immediately obtain \D I = ^A^3(012), implying triDi = -^d{d'^ - 1), 
since the dimension of the mixed symmetric subspace of (C'')®'^ is given by ^d{d^ — 1). 
The maximal averaged identification probability p^x^''(rf) for general d is then given by 

which is certainly less than the averaged maximum discrimination probability Pmax(c?) given in Eq.Q. 



IV. QUANTUM STATE IDENTIFICATION: CASE OF ARBITRARY iV 

The formulation proceeds in a similar way for arbitrary number (N) of copies of the reference states. We assume N 
copies of each reference state pa (a = 1, 2) are prepared in systems ai to a^. We denote the subsystem ai(8)a2<8' ■ ■'Sicin 
simply by a (a = 1, 2). The POVM now acts on the total system (g) 1 (g) 2 = (C"^)®(2Af+i) xhen the averaged success 
probability can be calculated in the same way as in the N = 1 case as follows: 

P^'^Hd) = I (tr[£;.p.(0)pi(l)®^p2(2)®^] 

-tr [Si(5Ar+i(01)5Ar(2) - 5Ar(l)5jV+l(02))] 

-tr 5^+1 (01)5^(2) - Sn{1)Sn+i{02) 



a=l,2 

1 1 

1 

2 



1 

< - 
- 2 



2djM^idN 
1 

4dN+idN 



id). 



(19) 



We note that it suffices to work in the subspace where systems 1 and 2 are both totally symmetric in order to 
determine non-zero eigenvalues and their multiplicities of the operator. 



D = 5Ar+i(Gl)5Ar(2) - Sn{1)Sn+i{02). 



(20) 



By the argument in the preceding section, we find the optimal POVM can be given by a projective measurement in 
this general case as well, Ei = Vd+ + 'Po and E2 — I — Ei. 
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A. Case of qubits (d — 2) 

In this case we find that the algebra of the angular momentum operators is useful. As stated above, we can work in 
the subspace where systems 1 and 2 are both totally symmetric, namely the total angular momentum of each system 
is The problem is then reduced to the recoupling of three angular momenta ^, ^ and ^. 

It should be noted that the support of the projector SN+i{Oa) (a = 1,2) is the subspace where the two angular 
momenta of systems and a are coupled to the total angular momentum + 5- Therefore, the projector 57v+i(0a) 
can be expressed in terms the angular momentum operators as 

SN+i{Oa) - (^2j{a) ■ s{0) + ^ + 1^ , (a = 1,2), (21) 

where s(0) is the spin operator of system and j{a) is the angular momentum operator of system a(= 1,2). We note 
2 V 2 



that j(a)' = f (f + !)■ 



Since 5Ar(a) is the identity in the subspace under consideration, we find 

2 



D 



N- 



-(i(l)-j(2)) •s(O). (22) 



Then can be readily calculated by use of the properties of the Pauli matrices and the commutation relations of 
angular momenta ji(l) and j{2). The result is given by 

where J — s{0) + j(l) + j(2) is the total angular momentum operator. From this we obtain the non-zero eigenvalues 
of D to be ±(5'-''^ with multiplicity 2 J + 1, where S'^"'"' is given by 



^N + lJ ' ' 2 2 
Thus we find the mean identification probability for qubit case to be 



1-(^V, iJ-il---,N-h. (24) 



pW(d^2)^- + - 1^ ^ y (2J+l)^/lT7Z±iV. (25) 

J^maxV ) 2 2(iV+l)(7V + 2) ^/ V^+ly 



1 1 

2+T^ 

1 1 

2 + 1^ 

1 1 





) for 


'rnax('^ 


= 2) 


(d 

max 


= 2) 


'rnax('^ 


= 2) 



2 , ^^VSc, 0.644, (26) 
2 , ^^{V2 + V5) ~ 0.703, (27) 



o ' on (^+4^3 + 3\/7) ~ 0.734 . (28) 
2 80 

The value of p[!^Jx{d = 2) in large N limit can be obtained by replacing the sum in H25|l with a continuous integration. 

(d = 2) i + ^ dxxv/l - x2 ^ I (iV ^ 00), (29) 
which is equal to the mean discrimination probability Pma.x{d = 2) as expected. 
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B. Case of arbitrary dimension d 

Now we study the case of general dimension, which is naturaUy more involved than the qubit case. Let us introduce 
the orthonormal base of the total system (C''^)8'(2JV+i) according to irreducible representations of the symmetric group 
S2N+1 and the unitary group U{d): 

I A, a, 6). (30) 

In this base, A represents an irreducible representation of S'2Ar+i, which is specified by a Young diagram. The 
expression A = [Ai, A2, • • •] means a Young diagram consisting of a set of rows with their lengths given by Ai, A2, • • •. 
The label a indexes vectors in a particular 52^+1 representation space and it runs from 1 to the dimension of the 
S'2Ar+i representation. The A also specifies irreducible representations of the unitary group U{d) and its vectors are 
indexed by b, which runs from 1 to mx{d), the multiphcity of representation A of S'2Ar+i on (C''')'^^^^^^^ 0|. We 
study the spectrum of D in the base of | A, a, 6). 

First we observe that D commutes with [/®(2^+i) fQj- arbitrary unitary U, since D involves only permutation 
operators among the 2N + 1 subsystems. By Shur's lemma, D should be proportional to the identity in a particular 
irreducible representation space of U{d). We can always choose the label a so that D is diagonal with respect to it. 
Then we have 

D\X,a,b)^A^JX,a,b}. (31) 

where is an eigenvalue of D, which is independent of b. 

We are interested only in non-zero eigenvalues of D, therefore, we work in the subspace, denoted by Vg, in which 
system 1 and 2 are both totally symmetric. The space Vg is the product of three U (d) irreducible representations 
[1] (g) [N] (X) [N]. Decomposing Vs into irreducible representations of U{d), we classify the base states in Vs into three 
groups. 

The first is the totally symmetric state: 

|[27V+1],6), 6=l,---,m[2jv+i](rf), (32) 

whose eigenvalues of D are obviously zero. The states in the second group belong to representations specified by 
Young diagrams of two rows [Ai, A2], where + 1 < Ai < 2N and Ai + A2 — 2N + 1. Since each of these U{d) 
representations appears twice in Vs, we distinguish the two by label a = ± as follows: 

I [Ai, A2],a,6), a = ±, b ^ 1, ■ ■ ■ ,m[Xi^x^]{d). (33) 

The remaining states are those whose Young diagram has three rows. We do not need the label a, since each of these 
representations occurs only once in Vg. 

|[Ai,A2,A3],5), b^l,---,m[x,MM]id). (34) 

Now we introduce the exchange operator T(12) between systems 1 and 2 by 

TV 

r(i2) = nr(i,2,). (35) 

1=1 

It is readily seen that T(12) anticommutes with D, implying that if | is an eigenstate of D, then T(12)| (j>) is also an 
eigenstate with the eigenvalue of opposite sign. On the other hand, the operation of T(12) does not change quantum 
numbers A or 6, since T(12) is just a permutation operator. Thus we conclude that the states with three rows have 
eigenvalue zero. And the label a = ± for the states with two rows can be chosen so that > and A^ = — A^. 
Therefore, tr|D| is given by 

2TV 

tr|7^| = 2 m[x,M]id)^+''''^- (36) 

Ai=TV+l 

We must still evaluate A'j^^''''^'. However, it should be noticed that this is independent of d and without loss of 
generality we can assume d = 2 in the calculation of a'^^^''^^'. This has already been calculated and given in Ea. H24() 
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in the preceding subsection. In the case of d = 2, the total angular momentum J is given by J = ^^2^^ = Xi — N — ^. 
We find 

(37) 



+ V + 1 



Thus we obtain the averaged success probability to be 



1 1 /A - /V\^ 

= 2 + 2^^ ^ 1 - [irTT) ' (^«) 

where m[x^^x^^{d), (A2 = 2A'' + 1 — Ai) is the multiplicity of the S'2jv+i irreducible representation [Ai, A2] in the total 
system 1 (g) 2, 

(Ai + d-l)!(A2 + ci-2)!(Ai-A2 + l) 
"^[^-^^l^'^) = (d-l)!(d-2)!(A: + l)!A2! ' ^''^ 

The optimal projective POVM can be expressed in terms of the states \X,a,b) as follows: 

2N 

E, = |[Ai,A2],+,6)([Ai,A2],+,6| + Po, 

Ai=Ar+l 

E2 = 1-Eu (40) 

where Pq is a projector onto any subspace spanned by zero-eigenvalue states of D, \ [2N +l],b) and | [Ai, A2, A3], 6). 

In the large A'' limit, we can verify that p^sL{d) approaches the averaged maximum discrimination probability 
Pmaxid) as follows: 



(41) 



V. CONCLUDING REMARKS 



Wc have discussed a natural generalization of pure state discrimination problem, which we call pure state identi- 
fication, where our task is to identify a given input state with one of the two reference states which arc presented 
as its N copies, but without classical information on them. The two reference states are assumed to be distributed 
in a unitary invariant way in the pure state space of dimension d. For arbitrary N and d, we have determined the 
averaged maximal success probability, which was shown to be achievable by a projective measurement. And it was 
explicitly shown that this mean identification probability approaches the mean discrimination probability in the large 
A'^ limit as expected. 

In this paper, we assumed that the input state is given as a single copy. If M copies of the input state are given, 

our task is to distinguish two states pf ^(0)pf ^(l)pf ^(2) and ^(0)pf ^(l)pf ^(2). By a similar argument to the 
one in the case of M = 1, the averaged maximum success probability can be shown to take the form: 

p'^^.fKd) = \+ . ■ V tr SN+MmSN{1) - SN{l)SN+Mm . (42) 

For general M,N, and d, however, it does not seem straightforward to evaluate this probability further. We would 
like to leave this problem to a future study. 
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